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Abstract 

We consider the asymptotic expansion of the sum 

p—an^ 

Sp{a;w) = Y,^^ 

n—1 

as a —>■ 0 in I arg a| < ^tt for arbitrary finite p > and w > 0. Our attention is 
concentrated mainly on the case when p and w are both even integers, where the ex¬ 
pansion consists of a finite algebraic expansion together with a sequence of increas¬ 
ingly subdominant exponential expansions. This exponentially small component 
produces a transformation for Sp{a; w) analogous to the well-known Poisson-Jacobi 
transformation for the sum with p = 2 and w = 0. Numerical results are given to 
illustrate the accuracy of the expansion obtained. 
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1. Introduction 

We consider the asymptotic expansion of the sum 

oo -anP 

Sp{a-,w) = Y,^^ ( 1 . 1 ) 

n=l 

as the parameter a —>• 0 in | arg a\ < ^tt, where p > 0 and, for convenience, w will be 
supposed throughout to be real and positive. When tc = 0, this sum is known as the 
Euler-Jacobi series and when a = 0 then S'p(0; w) reduces to the Riemann zeta function 
C(rc) (provided 3?(t(;) > 1). Consequently, the series in (1.1) can also be viewed as a 
smoothed Dirichlet series for C{w). 

The asymptotics of Sp{a]w) as a —)• 0-|- for p a rational fraction and w < 0 was 
considered by Ramanujan and is discussed in [1, Chapter 15]. It was shown that the 
expansion in this case consisted of an asymptotic sum involving the Riemann zeta 
function. A hypergeometric function approach for the Euler-Jacobi series when p is 
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a rational fraction and w = 0 has been discussed at length in the monograph [4], 
where great emphasis was placed on obtaining exponentially small expansions. The 
case re = 0 and arbitrary p > 0, thus generalising the work in [4], has been investigated 
in [6, §8.1] also for a —)■ 0+ using both a Mellin-Barnes integral approach and also a 
saddle point analysis of a Laplace-type integral representation. A similar Mellin-Barnes 
integral approach for this latter case has also been independently considered in [2]. The 
saddle point approach is instructive for understanding the appearance of exponentially 
small terms in the expansion as the parameter p increases. It was established in [6, 
§8.1] that an additional exponentially small contribution appears when p = 2 + Ak, 
k = 0,1, 2,... . Indeed, the appearance of the first exponentially small expansion as 
p passes through the ‘classical’ value p = 2 was demonstrated to be associated with a 
Stokes phenomenon [6, §8.1.7]. 

In the case p = 2, w = 0 the behaviour of the sum ^ 2 ( 0 ; 0) as a ^ 0 can be obtained 
from the classical Poisson-Jacobi transformation given by 

= = + (1.2) 

n=l n=l 

valid for all values of a in j arg a] < ^tt. This well-known transformation relates a sum 
of Gaussian exponentials involving the parameter a to a similar sum with parameter 
TT^/a. In the small-a limit, the convergence of the sum on the left-hand side becomes 
slow, whereas the sum on the right-hand side converges rapidly. Various proofs of (1.2) 
exist in the literature; see, for example, [6, p. 120], [8, p. 60] and [9, p. 124]. 

The dominant asymptotic expansion of Sp{a; rc) as a —)• 0 in j arg a] < \'k for general 
p > 0 and rc > 0 is relatively straightforward and is found to consist, in general, of a 
single term proportional to together with a series in ascending powers of a 

with coefficients involving the Riemann zeta function (the algebraic expansion). When 
0 < p < 1, the expansion is convergent and the result is exact; when p > 1 the 
expansion is asymptotic as a —)■ 0. The most interesting case arises when p and w are 
both even integers. The above-mentioned algebraic expansion then terminates after a 
finite number of terms, and it becomes essential for accurate estimation to also include 
a subdominant sequence of exponentially small expansions. This exponentially small 
component produces a transformation for Sp{a-,w) analogous to the Poisson-Jacobi 
transformation in (1.2), but valid as a —)• 0 in j arg a] < ^tt. This similarly involves a 
hnite sequence of series similar to (1.2) with a in the exponential replaced by an inverse 
power of a, but with each term decorated by an asymptotic series in ascending powers 
of 

The approach we employ in this paper is based on a Mellin-Barnes integral rep¬ 
resentation for Sp{a;w) and is similar to that described in [6, §8.1.4]. An algorithm 
for the determination of the coefficients in the exponentially small asymptotic series 
is described. The case p = 2 when w is an even integer has been recently discussed 
in [7], where the coefficients in the decorating asymptotic series can be given in closed 
form. An application of the series when p = 2, with w = 2 and w = 4, has arisen in 
the geological problem of thermochronometry in spherical geometry [10]. 


2. An expansion for Sp{a]w) as a —)■ 0 when w, p ^ 2,4,... 

We examine the expansion of the series Sp{a;w) defined in (1.1) as a —)■ 0 in the 
sector j arg a] < ^vr, where p > 0 and for simplicity in presentation we shall assume 
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throughout real values of w > 0. The case p = 1, w = 0 may be excluded from our 
consideration since the series in this case is summable as a geometric progression. 

Our starting point is the well-known Cahen-Mellin integral (see, for example, [6, 
§3.3.1]) 

rc+ooi 

z‘^e~^ = -—: / r(Q; — {z ^ 0, | arg z] < ivr), (2.1) 

ZTTl Jc-ooi 

where c < J?(q;) so that the integration path passes to the left of all the poles of r(a — s) 
situated at s = /c -|- a (A: = 0,1,2,...). Then, it follows that 


°° °° 'fl~^ r—c+ooi 


Sp{a;w) = = E 

n=l n=l 

r—c+ooi 

2>7ri J—c—ooi 


T{—s){an^yds 


27ri J —c—ooi 

T{—s)C,{w — ps)a‘^ds (| arg a| < iyr), (2.2) 


upon reversal of the order of summation and integration, which is justified when c > 
max{0, (rc — l)/p}, and evaluation of the inner sum in terms of the Riemann zeta 
function. 

The integrand in (2.2) possesses simple poles at s = /c (/c = 0,1, 2,...) and s = 
So = (ic — l)/p, except if sq = M (that is, w = pM J- 1), where M is a non-negative 
integer, when the pole at s = sq is double. The residue at the double pole is obtained 
by making use of the fact that C{s) ~ l/(s — 1) J- 7 in the neighbourhood of s = 1, 
where 7 is Euler’s constant, to find 

(—a+ (I 1 1 

\ 7 - - log a - y{M -M) I (M = 0,1, 2,...), 

M\ [ p p j 

where y(x) is the logarithmic derivative of the gamma function. The case when w and 
p are even positive integers requires a separate treatment which is discussed in Section 

3. 

2.1 The case 0 < p < 1 

We first consider the case 0 < p < 1. The integration path in (2.4) can be made to 
coincide with the imaginary s-axis together with a suitable indentation to lie to the left 
of the poles at s = 0 and s = sq (when 0 < in < 1). Then use of the functional relation 
for C(s) given by [9, p. 269] 

C(s) = 2^7r’^“^C(l — s)r(l — s) sin ^tts (2-3) 


shows that the integrand can be written as 


(27r)"'C(l- “ +ps) 


r(l — w + ps) sin ^T^ips — w) 


r(i + s) 


SmTTS 


(27r)P* 


With s = where i? —)■ 00 and is chosen so that the arc passes between the poles 
on the positive real axis, the logarithm of the dominant real part of the integrand is 
controlled by 

(p — l)Rcos6 log R + 0{R). 

When |6l| < ^tt, this last expression tends to —00 when 0 < p < 1. Consequently, the 
integration path can be bent back to enclose the poles of the integrand to yield the 
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convergent result 


Sp{a-,w) = Jp{a-,w) + - kp)a’^ 

k=o 


(0 <p < 1), 


where 


Jp{a;w) 


' -t(- — (wj^pM + l) 
P \ P J 

(—a)^ 1 1 

——{7 --log a + - V’(M +1)} {w=pM + l) 
\ Ml p p 


(2.4) 


and the prime on the sum over k denotes the omission of the term corresponding to 
k = M when w = pM + 1. 

When p = 1, the sum Si{a;w) is given by (2.4) but now the sum over k on the 
right-hand side converges when |a| < 27r; see [ 6 , §4.2.2] for details. If we let p = 1, 
w = 0 then use of the facts that C(1 ~ 2A:) = —B 2 k/{‘^k) and ({—2k) = 0, where B 2 k 
are even-order Bernoulli numbers, shows that 


5i(a;0) 


i _ 1 _i_ 1 „2A: 

a 2 ^a,^^( 2 fc)! 


(|a| < 27r), 


which correctly reduces to the trivial summation l/(e“ — 1) by application of [5, 
Eq. (24.2.1)]. 


2.2 The case p > 1 

When p > 1, the integration path in (2.2) cannot be bent back over the poles and we 
proceed in a similar manner to that described for the case re = 0 in [ 6 , §8.1.4]. Consider 
the integral taken round the rectangular contour with vertices at —c±iT, c'±iT, where 
o' > 0. The contribution from the upper and lower sides s = o" ± iT, —c < cr < c', 
vanishes as T —)■ 00 provided | arg a| < since from the behaviour 


T{a±it) = 0{t^ 2 e 2 ^*), ({a ± it) = log^ t) {t^oo), 


where for a and t real 


/i(fj) = 0 (cr > 1), 


^ — ^CJ (0 < fJ < 1), ^ — cr (cr < 0), 


T = l(0<cr<l), yl = 0 otherwise, 

the modulus of the integrand is controlled by 2 log Te~‘^^), with A = 

^TT — I arg a|. Displacement of the integration path to the right over a finite set of poles 
then yields (provided w and p are not even integers) 


N—l^ /_\fc 

Sp{a;w) = Jp{a;w) + '^'—r^C{w - kp)a^ + Rn, (2.5) 

fc=o 

where iV is a positive integer such that A > sq + f and the prime on the sum over k 
again denotes the omission of the term corresponding to k = M when w = pM + 1. 
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The remainder is given by 

rc-\-OQi 


pC+OOl 

Rn = 7 ;—: / T{—s)C,{w — ps)a^ds 

27r2 Jc—coi 


( 27 r)"' fc+ooi V(I - w + ps) sin ^ 7 r{ps - w) 

' C(1 “ + PS) —^- 


27ri 


r(i + s) 


smvrs 


(27r)P^ 


where c = N — ^ and the second expression follows from (2.3). Upon use of the result 
\C{x + iy)\ < ({x) when x > 1 , we obtain the bound 


\Rn\ < (27r)"' ^C{pN-p{so-^)) 

where cj) = arg a and 
F{t) = 


{2n)P 


N-i 


e-^^F{t) dt, 


( 2 . 6 ) 


r(i — w + ps) 


r(i + s) 


cosh ivrpt 


{s = N + 


cosh vrt 

On the integration path, F{t) is regular and satisfies F{t) = 0 (e“ 2 ^l*l) as t —)• ± 00 . 

2‘ 


Hence the integral in (2.6) is convergent and independent of |a| provided 101 < It 


then follows that 


Rn = 0{a^ 2 ) (a —)■ 0 in I arg a| < ivr). 


The expansion (2.5) is the dominant algebraic expansion associated with Sp{a;w) 
valid as a —)■ 0 in I arg a| < ^tt, provided w (> 0 ) and p are not even integers when the 
sum in (2.5) is finite. The same analysis can be applied to the case with non-positive 
w to yield the Berndt-Ramanujan result [1, Theorem 3.1, p. 306] 

5„(a;-u;) =^ ^ C{-w - kp) + 0 {a^-' 2 ) {w > 0 ) 

P \ P J 

as a —)■ 0 in I arg a| < ^tt. The reflection formula (2.3) can be employed to convert the 
argument of the zeta function to a positive form. 


3. The expansion of Sp{a;w) when w and p are even integers 

Throughout this section we let w and p be even positive integers, with w = 2m where 
m = 1,2,.... In this case sq = (2m. — l)/p, which cannot equal an integer and so no 
double pole can arise. More importantly, there is now only a finite set of poles of the 
integrand in (2.2) at s = sq and s = 0,1, 2,... , K, where K = [w/p\ , since the poles of 
r(— s) at s = ii' J- A: (/c = 1, 2,...) are cancelled by the trivial zeros of the zeta function 
(((s) at s = —2, —4,.... This has the consequence that the integrand is holomorphic 
in J?(s) > max{so,il"}, so that further displacement of the contour can produce no 
additional algebraic terms in the expansion of Sp{a;w). 

Thus, we hnd from (2.4) upon displacement of the integration path to the right over 
the poles of the integrand 

Spia; in) = -r f ^ _ pk) + (-)-(27r)-/i, (3.1) 

P \ P J ^ - 
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where 

II = ^ ^ r{-sK{w - ps)a^ds (3.2) 

and L denotes a path parallel to the imaginary axis with 3?(s) > (w/p) + S, with 6 
denoting an arbitrary positive quantity. This is easily seen to satisfy the requirement 
3ft(s) > max{so,R"o} necessary for the validity of (3.1). We now employ the functional 
relation for (((s) in (2.3) to convert the argument of the zeta function in (3.2) into one 
with real part greater than unity. The integral in (3.2) can then be written in the form 


II = 


2Tii 


C (1 — rc + ps) 


r(l — re + ps) sin ^vrps 


r(l + s) sinTTs 




where y = (27r)^a“^. 

In [ 6 , §8.1.4], the zeta function appearing in the above integrand was written as an 
infinite series. Here we follow a suggestion made in [2] and retain this function in the 
integrand; see also [4, § 8 ]. Making use of the expansion (see, for example, [ 6 , p. 368]) 


sm ^nps 
sinTTs 


AT-l 


= 2 ^ cos 7 r(|p— 2 r —l)s + 


r=0 


0 (p/2 even) 

1 (p/2 odd) ’ 


n = [Ip], 


where square brackets denote the nearest integer part^, we obtain 


N-l 

/l = ^ {J// + Jr} + 

r=0 


0 (p/2 even) 

J (p/2 odd) 


Here we have defined the integrals and J by 


1 


Jr = / C('i--w+Ps) 

2m Jl 


r(l-W+ps) 


and 


r(i + s) 




T(l-w + ps) 


r(i + s) 


(3.3) 

(3.4) 

(3.5) 


3.1 Asymptotie evaluation of and J 

The integrals and J have no poles in the half-plane 3ft(s) > (w/p) + d, so that we 
can displace the path L as far to the right as we please. On such a displaced path jsj 
is everywhere large. Let M denote an arbitrary positive integer. The ratio of gamma 
functions appearing in (3.4) and (3.5) may then be expanded by making use of the 
result (for p> 1) given in [ 6 , p. 53] 

^^r(lTy^ = ^(hK'^)-^!^Yl(-ycjT(KS + 'd-j)+pM(s)T(KS + 'd-M)y (3.6) 

where cq = 1, pm(s) = 0 ( 1 ) as jsj —)• oo in j arg sj < vr and 

K = p—l, h = p~P, ■d = ^ — w, A = (2 tt)2 k 2 ~'^p'^. (3.7) 

^The nearest integer part corresponds to [x] = N when x is in the interval [N — 1, -I- i]. Note 

that when p = 2, we have = 0 and the above expansion contains no information. 
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The coefficients Cj = Cj{w,p) (0 < j < 4) are listed in [6, pp. 46-48] where an algorithm 
for their determination is described; see Section 4 for details. 

Substitution of the expansion (3.6) into the integrals in (3.4) then produces 

Jt = ^ E C(1 - ^ + ps) r{KS + I? - i) (Xe^-^--)—ds + 7^± 

A ^^ f 

= ^ E j ^ Tiu)Ciqu + A,)(Xe^-^^)-“du + 7^± 

j=o ^ 

Here we have made the change of variable u —)■ Ks + 'd—j, with L' denoting the modified 
integration path, and have defined 

X := i>r := q ■=-, Xj := I +-{w + p{j - ^)) (3.8) 

K K K 

together with the remainders 

T^M,r = ^ij^ Pm(s)C(1 -w+ps) T{ks + I? - M) (3.9) 

We note that Xj > 0 for j > 0 when ze > 0 and p > 2. The above integrals appearing 
in may now be evaluated by means of (2.2), when we replace s by —s and allow 
the integration path (c — ooi,c + oof) to coincide with the path L', to yield the sum 
Sq{Xe^'"^^’’; Xj) as defined in (1.1). This evaluation is valid provided that the variable 
satisfies the convergence condition | arg(Xe^^*’^’')| < ^vr; that is 

—-— =F TT'ipr < At (0 < r < X — 1). 

K 

It is routine to verify that these conditions are met when | arg a| < ^vr. 

Thus we find 

A 

= Sq(Xe^^^^-,Xy + 7^E• (3-10) 

^ J=o 

Bounds for the remainders of the type have been considered in [6, p. 71, Lemma 

2.7]; see also [3, §10.1]. The integration path in (3.9) is such that J?(l — w + ps) > 1, 
so that we may employ the bound |C(a; + fy)| < C(x) for x > 1. A slight modification 
of Lemma 2.7 in [6, p. 71] then shows that 

K± = O (3,11) 

as a —)• 0 in the sector | arg a| < ^tt. 

An analogous procedure applied to J in (3.5) shows that 

4 

^ E (-)^c,X^-^Sq(X; X/) + 0(X^-^e-y (3.12) 

j=0 

as a —)■ 0 in I arg a| < ^tt. 

3.2 The expansion of Sp{a;w) 

The expansion of /l as a —)■ 0 in | arg a| < ^vr then follows from (3.1), (3.3), (3.10) and 
(3.12). We obtain the following theorem. 
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Theorem 1. Let m and M be positive integers. Then, when w = 2m and p is also an 
even positive integer, with K = \_w/p\ and N = [|p] {with square brackets denoting the 
nearest integer part), we have the expansion valid as a ^ 0 in \ arg a| < 


Sp{a-,w) = -r(" 
P \ 


1 — w 


D/P + ^ _ ph) a’^ + {-)^{27r)^h, 


k=0 


(3.13) 


with 


Af-l 

Il= Eria-,w,p) + Spp* EN{a;w,p), (3.14) 

r=0 


where Spp* is the Kronecker symbol with p* = 4iV + 2. The sums Er{a;w,p) are given 
by^ 


A 


M-l 


Er(a;w,p) = ^j:E Sg(Xe^^^^^;Xj) + RM,r 


(3.15) 


± j=o 


for 0 < r < N — 1, where X = K{h{2T:)P/a)^/^, ipr = {^p—2r — l)/K {0 < r < N — 1), 
Aj = 1 + {w + p{j — ^))/k, q = p/k and the parameters k, h, d and A are defined in 
(3.7). The leading coefficient cq = 1 and Cj = Cj{w,p) {j > 1) are discussed in Section 
4 . The sum E]\f{a]w,p) is also given by (3.13) when we put = 0 and omit the 
summation remainders RM,r satisfy the bound 

RM,r = (0 < r < A^). 


It is seen from (3.13), (3.14) and (3.15) that the sum Sp{a] w) has been expressed in 
terms of the sums Sq{Xe^'"^^’"Xj), which involve the reciprocal power of the asymptotic 
variable a scaling like Thus, as a —)■ 0 the argument X —)■ 00 . It is obvious from 

the definition in (1.1) (when q > 0) that 

Sq{z; Xj) ~ e~^ (z —)• 00 in I arg z\ < ^tt), 

so that the Er{a-, w,p) represent a series of exponentially small expansions of increasing 
subdominance in the small-a limit. In addition, the number of exponentially small 
expansions increases by one each time p increases by 4. By means of a saddle-point 
analysis in the case w = 0, this was demonstrated to correspond to a Stokes phe¬ 
nomenon when p was allowed to vary continuously through the values p = 2, 6,10 ,...; 
see [6, §§8.1.2, 8.1.7]. Finally, we remark that the exponents p and q are conjugate 
exponents [2], since 

1 1 

—i— — 1. 

p q 

When a is a real parameter, the expansion in Theorem 1 can be expressed in a 
different form by using (1.1) to represent the Sq{Xe^'^'''^^; Xj) as infinite sums. Then 
from (3.15) we obtain the following theorem: 

Theorem 2. Let w and p be even positive integers, N = [^p] and M be a positive 
integer. Then, the exponentially small expansions in (3.15) valid as a ^ 0-|- can be 
written in the form 

A °° 

Er(a;w,p) =— (0 < r < N - 1), (3.16) 


^The symbol signifies that the series with ± signs are to be added. 
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A °° 

EN{a-,w,p) = — (3.17) 

n=l 

The Tn^r {0 < r < N) have the asymptotic expansions 
M-l 

Tn,r = COs[Xn sill TTl/’r + Vr(j - I?)V’r] + 0{X~^), (3.18) 

j=0 

where Xn = Xn^^^ = K{h{2Trn)P/aY^'^, Yr = i^P — 2r — 1 )/n, Yn = 0 and the other 
quantities are as defined in Theorem 1. 

The result in (3.13) and (3.14) is the analogue of the Poisson-Jacobi transformation 
in (1.2) corresponding to re = 0, p = 2. In this latter case, iV = 0 and, from (3.14), 
II = Eo{a-, 0, 2). The ratio of gamma functions in (3.6) is replaced by the single gamma 
function r(s -|- ^) by the duplication formula for the gamma function, with the result 
that Co = 1, Cj = 0 (j > 1) and consequently Tn,N = 1 for all n > 1. Then, since K = 0 
and C(0) = —(3.13), (3.14) and (3.17) reduce to (1.2). The resulting expansion is 
valid for all values of the parameter a (not just a —)■ 0) satisfying | arg a| < ^tt. When 
w = 0, p = 2m, the expansions (3.16) and (3.17) reduce to those given in [7]; see 
Section 5. 

We observe that the n-dependence in the sums Er{a; w,p) from the factor n^~^X^ 
is given by Since p > 2 and tc > 0, this is seen to 

correspond to a negative power of n. 


4. The coefficients Cj 

We describe an algorithm for the computation of the coefficients Cj = Cj{w,p) that 
appear in the exponentially small expansions Ej.{a-,w,p) in (3.15) and (3.18). The 
expression for the ratio of two gamma functions in (3.6), with a=l — w for convenience, 
takes the form 

r(«+pg) ^ Cj , Pm{s) \ 

r(l J-s)r(«;s J-1?) 27r {l — KS — 'd)j {l — KS — 'd)M)^ 

where the parameters n, h, id and A are defined in (3.7) and {a)j = r(a + j)/T{a) 
is the Pochhammer symbol. If we introduce the scaled gamma function r*( 2 ;) = 
T{z)/{\/^ then we have 

r(/3s + 7 ) = T*{f3s + 7)(27r)ie-^^(/3s)^"+^-5 e(/Is; 7 ), 


where 

(/3s-b7- ^)log(l + ^) -7 • 
The above ratio of gamma functions may therefore be rewritten as 


e(fis;j) := exp 


R(s)G(s) 


M-l 


E 

j=0 


(1 — KS — 'd)j 


+ 


Pm{s) 

{1- KS- 'd)M 


(4.1) 



10 


R. B. Paris 


as |s| —)• oo in I arg s\ < vr, where 

e{ps-,a) 


Ris) = 


e(s; l)e(Ks; i?) ’ 


G{s) = 


r*(a +ps) 


r*(i + s ) t*{ks + ^)' 


We now let ^ ;= {ks) ^ and expand R{s) and G{s) for ^ 0 making use of the 

well-known expansion [ 6 , p. 71] 


k=0 


-k 


l^l ^ oo; I arg z\ < tt), 


139 


571 


2488320 ’ • • • 


where 7 ^ are the Stirling coefficients, with 

70 = 1) 7l = “12) 72 = 73 = 51840 ’ "14 

After some straightforward algebra we find that 

i?(,) = 1 + _ 1)| + 0(^2)^ = l + 0(^2)^ 

SO that upon equating coefficients of ^ in (4.1) we can obtain ci. The higher coefficients 
can be obtained by matching coefficients recursively with the aid of Mathematica to 
find [ 6 , p. 47] 

Co = 1, Cl = ■;^(2 — 5p + 2p^ — 12w + \2pw + I2w^), 


C 2 = „ (4 + 28p — 87p‘^ + 28p^ + 4p^ + 48i(; — 216pic + 24p^w + IMp^w 


I152p- 


—96t(;^ — 120pic^ -b 4:80p‘^w‘^ — 96u;^ + 480pu;^ -b 144u;^),.... (4.2) 

The rapidly increasing complexity of the coefficients with j > 3 prevents their presen¬ 
tation. However, this procedure is found to work well in specific cases when the various 
parameters have numerical values, where up to a maximum of 100 coefficients have 
been so calculated. In Table 1 we present the values^ of the coefficients Cj for 1 < j < 8 
in the specific examples considered in Section 5. 

When p = 2, use of the duplication formula shows that the ratio of gamma functions 
in (3.6) becomes 


r(i + s) 


j=0 


The coefficients cj in this case can be expressed in closed form as [ 6 , p. 53] 




ip = 2). 


(4.3) 


Finally, we mention that when w = 0 (corresponding to the Euler-Jacobi series) the 
coefficients Cj are listed for p > 2 and j < 8 in [ 6 , p. 374]. 


^In the tables we write the values as x{y) instead of a; x 10®'. 



Poisson-Jacobi-type transformation 


11 


Table 1; The coefficients Cj {1 < j < 8) for different p and w. 


J 

p = 4, w = 2 

p = 4, u; = 4 

p = 6, w = 2 

p = 6, w = 4 

1 

1.395833 (0) 

3.645833(0) 

1.472222 (0) 

3.305556(0) 

2 

3.495009(0) 

1.648980 (1) 

3.861497 (0) 

1.469946 (1) 

3 

1.230179 (1) 

9.075366 (1) 

1.380091 (1) 

8.081628 (1) 

4 

5.555372 (1) 

5.899040(2) 

6.207979 (1) 

5.260968(2) 

5 

3.060544 (2) 

4.424055 (3) 

3.387328 (2) 

3.949570(3) 

6 

1.990604(3) 

3.760330(4) 

2.188492 (3) 

3.358058(4) 

7 

1.493190 (4) 

3.572267(5) 

1.639364 (4) 

3.189927(5) 

8 

1.269216 (5) 

3.750863(6) 

1.396172 (5) 

3.348999(6) 


5. Numerical results and concluding remarks 


We present some examples of the expansion of Sp{a] w) given in Theorem 1 when p and 
w = 2m are even integers. For convenience in presentation, we extract the factor e~^ 
from the sum Sg{z; Xj) by writing 

OO -z(ni-l) 

Sg{z-, A,) = e-^S2{z; A,), A,) := ^ . 

n=l 


It follows that, when q > 0, Sq{z; Xj) = 0(1) as z —)• oo in | arg z\ < ^tt. 

Example 1. In the case p = 2, we have k = 1, q = 2, N = 0, K = m, p* = 2, i/jq = 0 
and X = tP' ja. The quantity bpp* = 1 so that the exponentially small component of 
^ 2 ( 0 ; 2m) consists of the single term Eo{a]2m,2). From (3.13), (3.14), (3.15) and (4.3) 
we therefore find 


S2{a]2m) - -F 


1 „/1- 2m 


a ™-2 


-E- 

k=0 


k! 


C(2m 


2h)a^ 


= i-r I r) 


2m-i 




1*^' (^)T2(^V«;2» + 2i) + 0(a«)} (5,1) 

as a —7- 0 in the sector | arg a| < ^tt. The expansion in this case has been given in an 
equivalent form in [7]. 

Example 2. When p = 4, we have k = 3, g = |,iV = l, iF = [|mj, p* = 6 and 


X = 3(^7r)^/^a , Xj = \{2m + Aj + 1). 

The quantity 5pp* = 0 so that there is the single exponentially small term Eq^g] 2m, 4) 
with i/iQ = |. Then we find the expansion as a —)■ 0 in | arg a| < given by 


54 ( 0 ; 2m) 


_ ^"^ ^Q(2m-l)/4 


K (_\k 

^^C(2m-4fe)a" 

k=o 
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= (-r 


9„\(4m—1)/6_ 1 . 


TT J 


M-1 




0(a^/3)| 

(5.2) 


where the coefficients cj = Cj{2m,A) can be obtained from (4.2) and [ 6 , p. 47] as 

Co = 1, Cl = ^(7 + 36m + 24m^), 

4o 

(385 + 4392m + 7104m2 + 3648m^ + 576m'‘), 

(39655 + 1191132m + 2970936m2 + 2666880m3 

+1080000m^ + 200448m^ + 13824m®),.... 


1 


C2 = 


C3 = 


4608 
1 

663552 


These coefficients are listed in Table 1 for 1 < j < 8 when m = 1 and m = 2. 

Example 3. When p = 6 , we have k = 5, g = |, = 1, RT = [^mj and 

X = 5(^71)®/® Aj = \{2m + 6 j + 2). 

In this case p* = 6 , so that 5pp* = 1 and there are now two exponentially small 

th 1 / 

^TT, we have the expansion 

K 


expansions £'o(a;2m,6), with i/’o = and Ei{a;2m,6). Then, as a —)■ 0 in | arg a| < 


SQ{a;2m) “ 


1/1 - 2 m 


6 


j(2m-i)/6 _ y- ^(2m - 6k)a’^ 


k=0 


k\ 


= {-r 


o^\(4m—1)/10 t- 2 ■ o 

4a V 


M-1 


TT J 


+i-y 


f IVl —. 

Sr 

S=0 (A 


(-)^9 


3^y4m-l)/10 


TT J 


3 

M-1 




)^' 


S6{Xe^5^^;Xj) + 


0(a^/®)l 


3=0 


0 (a^/®)l. 


(5.3) 


The first few coefficients cj = Cj(2m,6) are 


Co = 1, Cl = —(11 + 30m + 12m^), 
36 


1 


C2 = 


2592 

1 


C3 = 


1399680 


(517 + 3840m + 4116m2 + 1392m^ + 144m'^), 

(-22253 + 426550m + 8181720m2 + 5237640m3 
+1468800m^ + 185760m® + 8640m®) 


These coefficients are listed in Table 1 for 1 < j < 8 when m = 1 and m = 2. 

We show the results of numerical calculations to demonstrate the achievable accu¬ 
racy of the expansion in Theorem 1. We define the difference between Sp{a; w) and the 
finite algebraic expansion by 


Sp,w — ‘ 5 p^(fl) .— 5)5(0., rc) Tj 

p \ p 


{w-l)/p _ 


K 

E 

/c=0 


k\ 


C{vj—pk)a’^. (5.4) 
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Table 2: Values of the absolute error in the computation of Sp{a; w) defined in (5.4) using the expan¬ 
sions (5.1) and (5.2). The valne of the index jo corresponds to optimal truncation of the subdominant 
expansion Eo{a\w,p). 


a 

II 

to 

) " 

' to 

o 

o 

p = 2, w = 4 

Eq\ Jo 

1.00 

0.75 

0.50 

0.20 

0.10 

8.146(-06) 6.637(-09) 8 

2.031(-07) 8.089(-12) 11 

1.584(-10) 1.260(-17) 18 

5.774(-24) 1.542(-43) 47 

7.667(-46) 1.486(-86) 97 

6.252(-07) 3.642(-08) 6 

9.296(-09) 4.659(-ll) 9 

3.437(-12) 7.635(-17) 16 

2.189(-26) 9.830(-43) 45 

7.506(-49) 9.631(-86) 95 

a 

p = 4, w = 2 

SqI jo 

p = 4, w = 4 

j'o 

0.200 

0.100 

0.050 

0.010 

0.001 

3.473(-03) 1.329(-06) 7 

4.863(-04) 2.749(-08) 11 

2.737(-05) 2.156(-10) 14 

4.221(-09) 4.621(-17) 23 

1.064(-14) 1.033(-36) 53 

3.919(-04) 8.742(-06) 6 

4.805(-05) 4.879(-07) 8 

8.456(-06) 1.420(-09) 11 

7.982(-09) 3.041(-16) 21 

1.876(-16) 6.799(-36) 51 


In Table 2 we present the absolute error in the computation of Sp{a;w) for different 
values of the parameter a in the two cases p = 2 and p = 4, with w = 2 and w = 4 using 
the expansions given in (5.1) and (5.2). The first column in each entry displays the 
absolute value of Sp^w] that is, the accuracy achievable with just the algebraic expansion 
and no subdominant exponential terms. The second column shows the absolute error 
when the single optimally truncated exponential expansion EQ{a;w,p) (denoted by Eq 
in the table) is included. The optimal truncation index Jq, corresponding to truncation 
of the exponential expansion £'o(a;rc,p) at, or near, the least term in magnitude, is 
indicated in the final column. 

The situation when there is only a single subdominant exponentially small expan¬ 
sion present is straightforward: this sum is truncated at some suitable point thereby 
introducing a truncation error. If truncation is optimal, then the resulting error is 
exponentially more recessive than the parent exponential expansion. However, in the 
case of two, or more, exponential expansions of different degrees of subdominance (cor¬ 
responding to p > 6) the situation is not so obvious. It is not clear, without further 
investigation, how the error from the truncated leading exponential series compares 
with the contribution from the next series. 

We illustrate this by considering the case p = 6 and w = 2 given in (5.3). In 
Table 3 we present the absolute error in the computation of 56(a;2) as a function 
of the parameter a. We show, in order, the value of |56^2| and the absolute error in 
<^ 6,2 — £'o(tt;2,6) when the leading subdominant exponential expansion EQ{a]2,6) is 
optimally truncated at index Jq. The fourth column gives the absolute error when 
the hrst few terms of the second exponential expansion Ei{a;2,6) are included (for 
brevity in the table these exponential expansions are labelled Eq and Ei, and their 
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Table 3: Values of the absolute error in the computation of S'ela; 2) defined by (5.4) using the expan¬ 
sion (5.3). The value of the index jo corresponds to optimal truncation of the expansion Eo{a-, 2, 6). 


a 

‘^6,2 

1 x10"^ 

2.935(-02) 

5x10-2 

1.617(-03) 

1 x10-2 

9.512(-04) 

5x10-3 

1.292(-03) 

1x10-3 

1.604(-04) 

1 x10-^ 

9.894(-07) 

1x10-3 

6.209(-10) 


1^6,2 — Eq\ 

3.780(-05) 

3.037(-05) 

1.193(-07) 

1.099(-08) 

3.452(-ll) 

8.801(-17) 

1.522(-25) 


Jo 

o' 

1 

0 

Min £0 

6 

— 

9.422(-05) 

8 

1.200(-05) 

1.729(-05) 

12 

5.339(-08) 

1.228(-07) 

13 

8.713(-09) 

9.090(-09) 

19 

3.483(-12) 

3.757(-12) 

31 

2.230(-19) 

3.024(-19) 

51 

1.963(-30) 

1.964(-30) 


^i(j = 0) 


5.095(-05) 

1.191(-05) 

1.904(-07) 

2.148(-08) 

4.053(-ll) 

9.201(-17) 

1.564(-25) 


sum is denoted by Eq^i). The final two columns show the values of the least term 
(including prefactors) in EQ{a; 2, 6) at optimal truncation and the values of the leading 
term (j = 0) of the sub-sub dominant expansion Ei{a] 2, 6). 

A cursory inspection of Table 3 shows that for a ~ 0.1 the leading term of Ei (a; 2,6) 
is less than the minimum term of EqIu] 2,6) and consequently that inclusion of Ei (a; 2,6) 
cannot improve the accuracy. For a <0.01, the reverse is true: the leading terms of 
El (a; 2,6) are greater than the minimum term of Eo{a] 2,6) and their inclusion therefore 
increases the overall accuracy. However, it is clear that in both cases the final accuracy 
achievable is limited by the optimal truncation of the leading subdominant expansion 
£" 0 ( 0 ; 2,6). Further improvement in the accuracy would require a hyper asymptotic 
treatment in order to deal with the divergent tails of £ 0 ( 0 ; 2, 6) and £i(a;2,6). A 
possible hyper asymptotic scheme for the Euler-Jacobi series with p = 3 and w = 0 has 
been discussed in [4, §8]. 

Finally we remark that the asymptotics of the alternating version of (1.1) can be 
deduced from the result in Theorem 1 by making use of the identity 

00 p-anP 

-I 

n=l 
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